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ABSTRACT 

A  mathematical  continuum  theory  is  developed  to  describe  the  electromechanical 
behavior  of  dielectric  and  piezoelectric  solids  containing  imperfections.  The  macroscopic 
distortion  field  consists  of  recoverable  elasticity,  deviatoric  plasticity  arising  from  dislocation 
glide,  and  volumetric  deformation  from  vacancies.  A  connection  on  the  spatial  manifold  of 
deformed  lattice  vectors  describes  gradients  of  stretch  and  rotation  at  the  microscale  caused 
by  various  lattice  defects.  It  is  shown  that  parallel  transport  of  a  lattice  director  vector  with 
respect  to  this  connection  about  a  closed  loop  yields  a  discontinuity  with  contributions  from 
the  torsion  of  the  connection  (physically,  from  dislocations)  and  curvature  (physically,  from 
domain  walls  and  vacancy  gradients).  Classical  balance  laws  of  electrostatics  and  mass  and 
momentum  conservation  are  implemented.  A  free  energy  function  dependent  upon  lattice 
distortion,  polarization,  temperature,  and  defect  densities  is  suggested.  Thermodynamically 
consistent  kinetics  relations  for  dislocation  glide  and  vacancy  diffusion  are  then  derived,  with 
the  chemical  potential  for  the  latter  depending  upon  defect  density,  electric  potential, 
hydrostatic  pressure,  and  elastic  energy  density.  Vacancy  migration  and  mass  rearrangement 
at  the  free  surface  of  the  substance  are  also  considered  explicitly. 

Keywords:  dielectric,  piezoelectric,  diffusion,  vacancy,  dislocation,  disclination. 

1.  INTRODUCTION 

Dielectric  and  piezoelectric  solids  are  used  in  a  great  variety  of  solid  state  electronic 
devices  owing  to  their  electrical  and  mechanical  properties.  In  pure  dielectric  materials,  the 
mechanical  and  electrostatic  behaviors  are  generally  decoupled,  whereas  in  piezoelectric 
media  they  are  not.  In  the  latter,  mechanical  strains  may  induce  an  electric  field,  and  vice- 
versa.  Ferroelectric  media  undergo  a  change  in  crystal  structure  with  temperature,  and  may 
exhibit  dielectric  or  piezoelectric  effects  as  well  as  spontaneous  polarization  depending  upon 
crystallographic  phase.  Examples  of  ferroelectrics  encountered  in  the  semiconductor  industry 
include  the  crystalline  ceramics  PZT  (Lead  Zirconate  Titanate)  and  BST  (Barium  Strontium 
Titanate).  Defects  affect  the  reliability  of  many  ferroelectric  devices,  including  capacitors, 
high-frequency  sensors,  and  gate  dielectric  semiconductors,  notably  those  of  thin  film 
geometry  [1,2].  During  processing,  through  heat  treatments  or  addition  of  doping  chemicals, 
defect  densities  and  excess  electric  charges  may  be  controlled  [3].  Perhaps  the  most  naturally 
abundant  piezoelectric  material  is  crystalline  quartz  (Silicon  Dioxide).  Quartz  exhibits 
electromechanical  anisotropy,  with  historic  applications  in  pressure  transducers,  capacitors, 
and  gate  dielectrics  [4],  Quartz  is  also  a  fundamental  constituent  of  granite,  as  found  in  urban 
construction  applications  (e.g.,  buildings)  and  encountered  in  geological  settings  [5,  6]. 

Physical  models  are  needed  in  order  to  predict  how  existing  dielectric  and 
piezoelectric  materials  respond  to  new  environments,  to  design  engineering  systems 
comprised  of  such  materials,  and  to  tailor  new  materials  (e.g.  chemically  and/or  micro- 
structurally)  with  improved  combinations  of  properties.  Noteworthy  continuum  descriptions 
of  electromechanical  behavior  were  set  forth  in  the  mid-twentieth  century  [7-10].  More 
recently,  theories  based  on  thermodynamics  of  energy  minimization  have  been  developed  to 
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predict  domain  patterning  in  ferroelectric  media  [11],  and  vacancy  diffusion  and  surface 
rearrangement  in  dielectric  materials  [12,  13]. 

Historically,  defects  have  not  been  considered  explicitly  in  continuum  theories  of  such 
materials  [7-10],  even  though  their  presence  may  dominate  material  behavior  and  hence  the 
performance  of  the  engineering  system.  In  particular,  charged  vacancies  are  thought  to  affect 
leakage  current  characteristics  of  BST  thin  films  [3,  13].  Misoriented,  polarized  domains  are 
thought  to  influence  the  macroscopic  dielectric  constant  of  BST,  and  domain  walls  may  act  as 
a  source  or  sink  for  charges  and  point  and  line  defects.  Vacancy  migration  and  ion  tunneling 
are  important  phenomenon  in  quartz  resonators  [4].  Dislocations  have  been  found  in 
crystalline  quartz  [14]  and  their  motion  and  effects  on  properties  have  been  reported  [5,  6,  15- 
17].  Such  imperfections  affect  anisotropy,  piezoelectric  response,  and  mechanical  strength. 

Here  the  kinematic  description  of  defects  relies  on  methods  of  differential  geometry: 
the  purely  mathematical  basis  is  given  in  [18,  19],  while  recent  applications  to  mechanics  of 
solids  can  be  found  in  [20,  21],  In  the  present  work,  dislocations,  disclinations  (i.e.  rotational 
discontinuities)  and  point  defects  are  described  via  a  multiscale  continuum  framework 
following  [22],  The  disclination  concept,  used  previously  to  describe  subgrains  in  severely 
deforming  metallic  crystals  [23],  is  newly  applied  here  to  account  for  lattice  rotation  across 
domain  walls  in  polarizable  crystals.  The  description  relies  on  many  ideas  proposed  in  early 
field  theories  of  continuously  distributed  lattice  defects  [24-27]. 

In  the  present  theory,  the  mechanical  defonnation  is  considered  from  two 
perspectives.  From  the  macroscopic  perspective,  the  solid  deforms  according  to  the 
displacement  gradient  or  distortion,  with  contributions  from  elasticity,  plasticity,  and 
porosity.  Here,  distortions  are  assumed  small,  unlike  in  [22,  23],  as  fracture  is  thought  likely 
occur  in  the  relatively  brittle  materials  of  interest  prior  to  attainment  of  large  deformations. 
From  the  microscopic  perspective,  however,  deformations  may  be  large  in  the  vicinity  of 
defects,  including  significant  rotations  across  domain  walls.  A  linear  connection  is  defined, 
accounting  for  stretch  and  rotation  gradients  of  the  lattice  distinct  from  those  induced  by  the 
macroscopic  distortion.  The  total  discontinuity  or  burgers  vector  in  the  lattice  arising  from 
this  connection  is  shown  to  consist  of  contributions  from  dislocations,  disclinations,  and 
vacancy  gradients.  Bulk  and  surface  fluxes  of  vacancies  compatible  with  mass  conservation 
principles  are  introduced  [12,  13,  28].  Momentum  and  energy  conservation  laws  are  defined, 
along  with  standard  relations  from  electrostatics  [7].  Thermodynamic  dependencies  and 
kinetic  relations  are  then  postulated,  describing,  respectively,  the  state  of  the  substance  and 
the  motion  of  defects.  A  proof  is  given  that  such  relations  are  consistent  with  laws  of 
thermodynamics,  inspired  by  procedures  of  [29-31].  Finally,  a  specific  functional  form  is 
suggested  for  the  free  energy  density  of  dielectric  and  piezoelectric  bodies,  leading  to  a 
particular  kinetic  law  for  charged  vacancy  diffusion. 

2.  KINEMATICS  AND  BALANCE  LAWS 

The  mechanical  and  electrical  displacement  vectors,  ua  and  D" ,  are  defined  by 

ua=xa-xa0,  (1) 

Da  -  s0E“  +  P“ ,  (2) 

where  xa  and  x‘‘  are  deformed  and  initial  material  coordinates,  Pa  is  the  polarization 
defined  only  within  the  solid  domain,  s0  is  the  permittivity  of  free  space,  and  Ea  is  the 
electric  field  satisfying 

£„=-*«.  (3) 
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with  (j)  the  electric  potential  [7].  Indices  in  Roman  font  span  three  spatial  coordinates,  and 
subscripted  commas  denote  spatial  derivatives  with  respect  to  these  coordinates.  Polarization 
implies  a  relative  displacement  among  charged  ions  comprising  the  primitive  unit  cell  of  the 
lattice.  The  displacement  gradient  is  decomposed  additively  as 

Ua,b  =  Tab  +  Yab  +  Y'ab  ’  (4) 

with  ylh  the  recoverable  lattice  distortion,  yPab  the  plastic  distortion  from  dislocation  glide, 
and  yvab  expansion  or  contraction  from  point  or  volumetric  defects,  the  latter  two  satisfying 

(/)>  °.  (/)•„={  x/ysi,  (5) 

where  x  is  the  volume  fraction  of  defects  and  d"h  is  Kronecker's  delta.  Plastic  flow  is 
dictated  by  [23,  32,  33] 

fab=Tjv{k)sa)mlk)  ’  (6) 

k 

where  v  is  the  slip  rate  on  glide  system  k,  with  slip  direction  sa  and  slip  plane  normal  mb . 
Note  that  sa  and  mb  are  orthogonal  by  definition.  Defects,  polarization,  and  slip  plane 
geometry  are  illustrated  in  Fig.  1. 


Fig.  1.  Material  element  with  point  and  line  defects,  polarized  domains,  and  slip  plane. 

Diffusion  of  point  defects  is  specialized  here  to  vacancies.  Let 

>  (7) 

where  E,  is  the  number  of  vacancies  per  unit  volume  and  a  is  a  positive  constant.  Then 

vacancy  conservation,  in  the  absence  of  generation  or  annihilation,  requires 

i=-o:a,  (8) 

where  Q‘  is  the  bulk  flux  of  vacancies.  Surface  fluxes  are  also  admitted,  described  by  the 
vector  qa ,  where  Greek  indices  span  two  surface  coordinates.  Let  a  denote  the  nonnal 
velocity  of  the  surface,  which  may  contract  as  vacancies  are  released.  Denote  by  M  the  total 
mass  of  the  system,  with  actual  and  referential  (undeformed  and  defect-free)  mass  densities 
p  and  Pq  .  Conservation  at  surface  s  of  the  system  requires  a  null  sum  of  rates  of  change  in 
mass  from  surface  expansion/contraction  and  from  bulk  and  surface  fluxes  of  vacancies: 

M  /  p0  =  \a(l-x~rEaa  )ds  +  a\{Q“na  +qaa)ds  =  0,  (9) 

leading  to  [12,  13] 
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i=-a(Q‘na+ql)l(l-x-r«)- 


(10) 


Bulk  and  surface  coordinates  and  fluxes  are  illustrated  in  Fig.  2. 


Fig.  2.  Defect  fluxes  and  coordinate  systems. 


The  lattice  deformation  Aab  is  introduced  as 


<=#:+rEab+/ 


.b  9 


(11) 


and  includes  effects  of  elasticity  and  stress-free  volumetric  expansion,  but  not  dislocation 
glide,  as  the  latter  presumably  does  not  alter  the  crystal  structure  once  dislocations  have 
convected  through  the  material  [22,  23],  Lattice  director  vectors  da  (see  Fig.  1)  are 


deformed  via  (11)  to  the  vectors  d bAa  ,  and  the  corresponding  metric  for  their  change  in 
length  is  given  by  the  inner  product 

(12) 


Cab=A^dc-A-lddd=A: 


cScdA~ 


where  the  da  constitute  an  initially  orthogonal  triad  attached  to  each  unit  cell  of  the  lattice, 

and  do  not  necessarily  correspond  to  interatomic  bond  vectors.  Spatial  gradients  on  the 
manifold  of  lattice  directors  are  described  by  the  covariant  derivative  [18,  19,  22] 


VA=d  a,b-rbcadc, 

where  Christoffel  symbols  of  the  connection 


r  a  -  Aa  A-,d  +  +  r  Aa 

1  cb  ~  ,b,c  +  ^cb  +  1  c°b  • 


(13) 

(14) 


Above,  E'cl  and  T c  are  additional  degrees  of  freedom  accounting  for  rotation  and  stretch 

gradients,  respectively,  due  to  surface  and  point  defects.  Covariant  versions  of  these 
quantities  are 


J  cba 


Z7-dr  =  « 

“  cb  ^ da  “c[6a]  ’ 


rcsbdcda  =rcc(ba), 


(15) 


where  bracketed  indices  are  anti-symmetric  and  indices  in  parentheses  are  symmetric.  The 
nonvanishing  components  of  the  torsion  T  and  curvature  R  tensors  of  the  connection  (14)  are 
then  found  as  [19,  22,  27]: 


*K=  2^“1M+2^]+2r[c^, 

(16) 

r^\cd\~^ e\ba\  5 

(17) 

R{ab)[cd]  =  ^[c^d]^(ab)  +  WeC{ab)  • 

(18) 

The  microstretch  renders  the  connection  (14)  non-metric,  since  V(.C.(/)  =  -2YcCab. 
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Definition  2.1  Dislocation  tensor, 
form  as 


The  dislocation  density  tensor  is  defined  in  component 


/-> — ab  bed  rjn..a 

2a  =  e  Tdc 


(19) 


Remark.  When  the  microrotation  and  microstretch  vanish,  aab  in  (19)  agrees  with  the  form 
given  in  [33]  in  the  absence  of  volumetric  defects: 


—ab  ~ 

a  -  s 


bed 


AaA~\  „  = 


bed  a  a  A-le 
6  A. e.d A  .< 


bed  .Ea 

-e  y  x 


,d  9 


(20) 


and  vanishes  when  the  elastic  distortion  is  integrable,  i.e.  when  yEab  =  uEah ,  where  it '  is  an 
elastic  displacement.  In  (19)  and  (20),  sabc  denotes  the  permutation  operator. 


Definition  2.2  Disclination  tensor.  The  rank  two  disclination  density  is  defined  as 

4  9sd  =s8ba£dceRabce.  (21) 

Remark.  The  disclination  density  vanishes  when  the  microrotation  Ecba  =  0  . 


Definition  2.3  Point  defect  vector.  This  quantity  is  defined  by 

»’  =  c”"V"'y„,w  =  zs-'Vj  rd] + r, .  (22) 

Remark.  The  point  defect  vector  vanishes  when  the  microstretch  T c  =  0  . 


Definition  2.4  Total  burgers  vector.  The  total  burgers  vector  over  oriented  lattice  plane  a 
from  defect  fields  represented  in  (14)  is  given  by 

ba  =\(aad  +  C~Xaf£fgbe  sdxh  +  3d xa ) ndda  .  (23) 

a 

Remark.  In  the  absence  of  microstretch  and  microrotation,  alb  can  be  interpreted  as  the 
number  of  parallel  dislocation  lines  per  unit  area  with  constant  burgers  vector  b 
perpendicular  to  oriented  differential  area  element  nda  . 

Proposition  2.5  The  total  burgers  vector  is  equivalent  to  the  following  discontinuity  induced 
by  parallel  transport  of  a  lattice  vector  about  an  infinitesimal  closed  circuit  c: 

ba  =  \§(Tbc  - Kd*  ) dxb  a dxc  .  (24) 

c 

Proof.  The  contribution  from  torsion  and  curvature  to  the  closure  failure  upon  parallel 
transport  with  respect  to  a  connection  of  the  form  (13)  is  derived  in  Eqs.  (2.19)  and  (4.1)  and 
accompanying  discussion  of  [19],  neglecting  terms  of  third  and  higher  order  in  dx .  Then, 
from  the  skew-symmetry  of  the  differential  two-fonn  dxh  a  dx  =  £dbcndda  ,  and  substituting 
definitions  (19)-(22)  into  (24),  expression  (23)  is  recovered. 

The  remaining  conservation  laws  are  now  introduced  and  are  used  in  the 
thennodynamic  description  that  follows.  Bulk  mass  conservation  is  described  by 

P  +  PK=  0-  (25) 


5 


Linear  and  angular  momentum  conservation  are  written 

aabb  =  0 ,  (7[ab]=0,  (26) 

and  the  mechanical  traction  vector  ta  =  crahnh  is  introduced.  Maxwell's  equation  of 
electrostatics  is 

P  =  Da,a,  (27) 

with  the  charge  density  related  to  volumetric  defects  (i.e.,  charged  vacancies)  as  follows: 


p  =  ez%. 


(28) 


where  e  is  the  charge  of  an  electron  and  z  is  the  valence  contribution  of  each  defect. 


3.  THERMODYNAMICS  AND  KINETICS 

Two  energy  balances  are  introduced.  The  first  is  the  local  balance  of  electrostatic 
energy  [7,  30]: 

~(e„E‘  +  P%t  =  fp^  D‘Ea  =  fp .  (29) 

The  second  is  the  global  balance  of  the  rate  of  external  work,  internal  energy,  and 
electrostatic  field  energy  [13,  30]: 

W  =  f  pUdv  +  —  [  —EaEadv  +  —  [  (7<f)ds  ,  (30) 

J  dtJ  2  ‘  dt J 

where  U  is  the  internal  energy  per  unit  mass,  &  is  the  surface  charge  density,  and  the  rate  of 
external  work  W  is  given  by 

W  =  J  <j(/)ds  +  J  tauads  -  J  juQanads  ,  (31) 

with  ju  the  chemical  potential  for  vacancy  diffusion.  The  Helmholtz  free  energy  <// ,  specific 
entropy  rj ,  and  absolute  temperature  0  are  related  by 


II 

c; 

i 

-cs 

<$> 

(32) 

Neglecting  heat  conduction,  the  global  entropy  inequality  is  written 

^  Opfjdv  >  0  ^  ^  p{ll  -y/  -di{^dv>0  , 

(33) 

and  can  be  expressed  compactly  as 

o  =  w-n>  o, 

(34) 

where  the  rate  of  system  energy 

&  =  J P (w  +  0*! ) dv  +  ^  ^  ^  2  dv  +  J &</>ds  . 

(35) 

The  Helmholtz  free  energy  is  assumed  to  exhibit  the  following  functional 
dependencies: 

V'=V'(r‘d,,P‘,  f.M),  (36) 
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with  the  generic  scalar-valued  function  h  =  h(aab ,0ab )  accounting  for  energetic 
contributions  from  defects  apart  from  E,  • 

Proposition  3.1  A  set  of  bulk  thermodynamic  and  kinetic  relations  simultaneously  satisfying 
(29),  (34),  and  constitutive  assumption  (36)  and  neglecting  electrostatic  terms  of  second  order 
is 


°"b=pdvE- 

byEab 

(37) 

p  _  dip 
a  p  dPa  ’ 

(38) 

dw 

1  = - > 

do 

(39) 

Qa  =  -dah ph  \dab\>  0, 

(40) 

dip 

P  =  P—  +  ez<j)  +  ap, 
d£, 

(41) 

v0  rw/cr0  "sgn(rw),  v0>0, 

(42) 

Ytmv{k)>p^-h. 
k  dh 

(43) 

Above,  dab  is  the  diffusivity,  T<k)  =  (Jahs{k)m{k)  is  the  resolved  shear  stress  on  slip  system  k, 
v0  is  a  reference  shearing  rate,  n  denotes  rate  hardening,  and  a0  is  the  slip  resistance.  The 
notation  sgn  jc  =  x  /  |jc|  . 


Proof.  Expanding  the  internal  energy  rate  using  (36)  gives 

p{v  +  Ori)  =  P 


dyt 


ab 


dPa 


dh 


dw 

—!—  +  ri 
ySO  J 


e. 


(44) 


from  which  assuming  that  Q  -  0  for  purely  thermal  processes  gives  (39).  Expanding  the 
electric  field  energy  by  using  Eq.  (29), 


E-EJv  ]  =  -J  E»P-dv  + 1  +  J  u-I-^dv -  J  r\utds ,  (45) 


where  the  divergence  of  the  Maxwell  stress  [9,  30], 

-A  =  ~tP  =  (E.D>  ~ /  2)4  • 


(46) 


as  well  as  the  value  of  the  stress  itself,  Sah ,  are  both  set  to  zero  henceforth  by  the  assumption 
of  small  electrostatic  terms.  The  net  dissipation  of  (34)  then  becomes 


0  =  J  tauads  -  J  paQadv  -jp  yEabdv 

^  t  ab 

\{E.-P^)rdvf{p^+ez4,-p^i;dv-\p^hdv. 


(47) 
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Expanding  the  stress  power  using  the  divergence  theorem,  (4),  and  (5), 

I  fu,ds  =  J<t£  +  <x“»  (ft  +  y'  +  a&t  /3)rfv  =  J (CT“V*  +  <r-v:4  -  pa|)rfv ,  (48) 

0 

where  the  deviatoric  stress  dab  =  <jab  +  pSab  and  the  hydrostatic  pressure  -3 p  =  cr“a .  Then 

0  =  J  VabrPabdv  -  J  M,aQ“dv  ~  j  p  ^ hdv  +  ”  ' • 


^  -  P^yEahdV-\[Ea-  P^)padv~\yp^  +  ez(l)- ^  +  pa)^ 


dip 


dy/ 


(49) 


The  final  three  volume  integrals  on  the  right  of  (49)  vanish  when  (37),  (38),  and  (41)  apply. 
From  (40),  dissipation  is  always  non-negative  from  vacancy  flux: 


-paQa=padabpb>  0.  (50) 

This  leaves  the  local  requirement 

®"hyFab  -  p  -  o  •  (si) 

oh 

From  (42),  plastic  work  is  always  non-negative: 

GiJY?j  =  X r a)vo  |r W  /  °o I"  s§n {T<k) )  -  0  ’  (52) 

k 

meaning  that  the  energy  storage  rate  from  generation,  annihilation,  and/or  interaction  of 
lattice  defects,  h ,  cannot  exceed  the  rate  of  dissipation  from  dislocation  glide,  as  required  by 
(43). 


Corollary  3.2  Thermodynamically  consistent  boundary  conditions  on  bulk  vacancy  flux,  and 
an  admissible  kinetic  equation  for  surface  flux,  compatible  with  Proposition  3.1,  are 


assuming  that 


Q\  =  fi 


a 


{py/  +  (j)p)l(\-X-YEaa)  +  p 
{p¥  +  <l>P)l{\-X-YEaa)\ 


p>0, 

>0. 


(53) 

(54) 


(pip  +  </>p)>0,  AaP\{p\i/  +  fp)l{\-x-YE“ )]^  >0, 


(55) 


and  subject  to  the  boundary  restrictions  of  constant  (in  time)  charge  and  electric  potential  on 
s,  and  traction-free  or  quasi-static  surface  s. 


Proof.  The  system’s  total  energy  rate  with  a  moving  boundary  is  written  using  (35)  and  (45): 

f2  =  ^ p{y  +  Or/^dv- ^ EiP'dv  +  J (j>pdv  d</>ds  +  J a(py/  +  ez<f>^ ] ds  ,  (56) 

assuming  Maxwell  stress  Iab  =  0.  This  is  expanded  using  (36)  and  (48)  as 


Q  =  | tauads - J dabyPahdv  +  J p^-hdv+jyp^^  +  ezfi  +  ap  |Jv  +  | ai^py/  +  ez(j)f}ds  .  (57) 
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For  a  traction-free  or  quasi-static  surface,  taua  =  0 ,  and  with  dissipative  plastic  work  and 
plastic  variables,  (5 1)  applies.  Then  using  the  divergence  theorem, 


Q  =  -J juQanads  +  J p  aQadv  +  j  a(pi//  +  ez(f>^ds  .  (58) 

Since  jU  aQa  =  -dab jua/uh  <  0,  the  following  requirement  arises  in  order  to  satisfy  (34)  in  the 
absence  of  external  working: 

d(pi//  +  </>p)-juQana<0.  (59) 

Substituting  from  mass  conservation  requirement  (10)  gives 

-Qana[a{p\//  +  (j)p)l[l-x-YEaa)  + p\  +  qaa(pW  +  <l>p)l(\-X-rEaa  )^°-  (60) 

Then,  from  (53), 

-QanJ\jx{pi//  +  (j)p)l[\- x~yEaa)  +  p\  =  -p\a{py/  +  (/)p)l[\- x~rEa)  +  m\  ^0.  (61) 

Finally,  from  (54)  and  (55),  the  following  inequality  remains: 

qaApv/+(l>p)l(x-x-rEaa)  = 

-  AaP  \[PY  +  (/>p)l{  l-Z-  yEaa )]  Pa  [PY  +  $>)■ /  (l  -  -  rEa  )  ^  0. 

Remark.  Since  0  <  %  <  1 ,  (55)  and  (62)  suggest  that  the  sum  of  the  free  energy  per  unit 
volume  and  electrostatic  energy  per  unit  volume  should  be  non-negative  and  convex  when 
elastic  volume  changes  are  small. 

4.  CONSTITUTIVE  FUNCTIONS  FOR  PIEZOELECTRICS  AND  DIELECTRICS 

The  following  more  specific  form  of  (36)  is  proposed  for  single  crystalline 
electromechanical  solids,  in  the  absence  of  phase  changes  or  spontaneous  polarization: 

P^  =  (i/2)c^(i-^)r(VM+(1/2K^  +  ^r,^  +  ^.«)+A.  («) 

where  C  and  A  are,  respectively,  the  anisotropic  elastic  moduli  and  inverse  dielectric 
susceptibility,  and  A  are  piezoelectric  constants.  For  a  dielectric  material  such  as  BST  in  the 
paraelectric  state  [3,  13],  Aab  =  0 ,  and  (2),  (37),  and  (38)  lead  to 
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(67) 

(68) 
(69) 


where  I7ahc  =  A  lad Ahcd  .  The  scalar  function  cp  is  assumed  to  follow  the  universal  relation  for 
the  chemical  potential  of  an  ideal  mixture  [13,  34]: 

<PN A  =  Nt(G0  (0)  +  NAkBOv  In  v) ,  (70) 

where  NA  is  Avagadro’s  number,  kB  is  Boltzmann's  constant,  Nr  =1/  a  is  the  number  of 
atomic  sites  per  unit  volume,  v  =  % /  Nr  is  the  mole  fraction  of  vacancies,  and  G0  (6  )  is  the 

bulk  Gibbs  free  energy  of  the  unstressed,  defect-free  crystal.  Note  that  for  constant 
temperature  applications,  each  atom  (or  missing  atom)  occupies  a  fixed  reference  volume, 
such  that  Nt  =  1  la  and  v  =  x  ■  From  (41),  the  chemical  potential  for  vacancy  diffusion  is 
then 

n  =  kte(  1  +  In  X) + ezj  -  a  [(1  -  X)  /  3 + (1  /  2)  C ]  •  (71) 

The  kinetic  equation  for  bulk  flux  follows  from  (40)  as 

Q‘  =  d*  \  <*Eb  -k-A-Xt+  ^0(c  frU„  -  ) + “C*  A)A).»  1  •  <72> 

X  -7 

Consider  now  a  one-dimensional  reduction  of  (72),  where  uniaxial  strain  conditions  hold  for 
the  x-direction,  i.e.  a  =  C(l-/)f£  ,  where  C  is  an  elastic  constant: 


kB0  dx  a(\-x) 


Q/d=ezE — ^ - ^  + 

X  dx 


C(l  +  3r£)- 


-c  rE  — 

dx 


This  can  also  be  written  as  follows,  neglecting  terms  of  second  order  in  elastic  strain  and 
porosity: 


Q!d=- 


2aCyE 


kBe  dx 


-[ez]^j—[a(.1- 2r)(1+3r£)]- 


where  the  first  tenn  on  the  right  causes  vacancies  to  move  to  areas  of  lower  concentration,  the 
second  term  causes  charged  vacancies  to  move  to  areas  of  lower  electric  potential,  and  the 
third  tenn  accounts  for  the  elastic  strain  energy  density,  causing  vacancies  to  move  to  areas  of 
lower  hydrostatic  pressure.  While  the  final  term  vanishes  under  static  stress  equilibrium  in 
the  one  dimensional  case  (i.e.,  the  first  of  (26)  gives  dp  /  dx  =  0),  analogous  terms  containing 
stress  or  pressure  gradients  do  not  generally  vanish  in  the  three  dimensional  Eq.  (72). 


5.  CONCLUSIONS 

A  continuum  framework  for  the  electromechanical  behavior  of  dielectric  and 
piezoelectric  solids  has  been  developed.  The  kinematic  description  accounts  for  the 
contributions  of  dislocations,  disclinations  (i.e.,  domain  walls),  and  point  defects  to  covariant 
derivatives  and  discontinuities  in  the  field  of  lattice  director  vectors  on  their  spatial  manifold. 
Thermodynamically  admissible  relationships  for  recoverable  elasticity,  polarization, 
dislocation  glide,  and  bulk  and  surface  diffusion  of  charged  vacancies  have  been  derived. 
Specific  constitutive  functions  have  been  postulated  for  dielectric  and  piezoelectric  single 
crystals,  with  particular  emphasis  given  to  the  form  of  the  governing  equation  for  bulk 
diffusion  kinetics. 
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